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^ ■ Abstract. We study in detail the effect of quasicondensation. We show that this 

effect is strictly related to dimensionality of the system. It is present in one dimensional 
systems independently of interactions - exists in repulsive, attractive or in non- 

\^ i interacting Bose gas in some range of temperatures below characteristic temperature of 

the quantum degeneracy. Based on this observation we analyze the quasicondensation 
in terms of a ratio of the two largest eigenvalues of the single particle density matrix 

OO , for the ideal gas. We show that in the thermodynamic limit in higher dimensions the 

second largest eigenvalue vanishes (as compared to the first one) with total number of 

rij ■ particles as ~ N^'^ whereas goes to zero only logarithmically in one dimension. We 

i^ , also study the effect of quasicondensation for various geometries of the system: from 

quasi-lD elongated one, through spherically symmetric 3D case to quasi-2D pancake- 
like geometry. 

C^ . PACS numbers: 67.85.Bc, 03.75.Hh, 05.30.Jp 
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1. Introduction 

Properties of tlie interacting, harmonically trapped, ultracold gas are much more 
interesting in one dimension than in two and three dimensions. In 2D and 3D a 
phase transition occurs whereas in ID does not (surprisingly in quasi- ID systems two- 
step condensation is possible jM])- Moreover, in symmetric 3D systems it was shown 
experimentally [22| [TT] that nearly all the way up to critical temperature the phase 
of the cloud is spatially uniform, equivalently: a coherence length of the system is 
equal to its size. On the contrary, Petrov et al. have shown that the observation 
of the quasicondensate (condensate with fluctuating phase) is possible in ID [19], 2D 
|18| and very elongated 3D [20] repulsive Bose gas. Those predictions were confirmed 
experimentally [9l [211 [12] • Recently density and phase properties of elongated systems 
were investigated in a number of experiments: [TOl [23l [3l [I5] . In spite of many theoretical 
attempts [2], [H [TJ [8] the theory of quasicondesates is still not as mature as the theory of 
condensates. In this paper we shed a new light on the quasicondensation phenomenon. 
In [6] we stressed that shortening of the coherence length is not only the property of a 
repulsive gas but also of an attractive one. There exists a direct connection between the 
quasicondensation and the spectrum of a one-body density matrix [j. More precisely, the 
quasicondensation occurs when occupation of more than one eigenmode is comparable 
with the total population ^7\. In some sense it is similar to the "fragmented" spinor 
condensate [1^ . In Section [2] we show that properties of one dimensional Bose gas 
are similar regardless the sign of the interaction. In Section [3] we explore coherence 
properties of the ideal gas in ID, 2D, 3D and with arbitrary ratio of the trap frequencies. 

2. Interacting gas 

Firstly, we study a one dimensional, weakly interacting Bose gas confined in a harmonic 
trap. Thus, our Hamiltonian of the one dimensional Bose gas has a form: 

H = ^^(x) I ^ + -mwlx^ j i!{x)dx + 

+ - l¥{x)¥{x)ii{x)ii{x). (1) 

The Hamiltonian is a sum of the single particle oscillator energy with mass m and 
angular frequency uj and a conventional contact interaction with the coupling constant 
g. Throughout this paper we use the oscillator units of position, energy and temperature. 



^Wo and -^ respectively. Hence a dimensionless coupling g is in units of 



mwQ ' k 



Almost all results are calculated for 1000 atoms. 

At the beginning we study a one dimensional, repulsive, weakly interacting Bose 
gas confined in a harmonic trap. Our results are for the canonical statistical ensemble, 

\ In the description of partially coherent light spatial coherence modes are used also as eigenfunctions 
of the first order correlation function of the light field. See for instance B. Saleh Photoelectron Statistics: 
With Applications to Spectroscopy and Optical Communication (Springer- Verlag, New York, 1978) 
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Figure 1. (Color online) Correlation length, width of g\ and size of the whole cloud. 
All lengths are defined as full width at half maximum. Results for (a) repulsive and 
(b) attractive gas. 1000 atoms considered. 



thus the temperature is a control parameter. In previous works we showed that in a 
wide range of temperatures such a system can be efficiently described using the so called 
classical field approximation O [5] . To analyze the coherence properties of the system 
we use the first-order correlation function defined as: 

(^*(-x)^(x)) 



gi{-x,x) 



(l^(^) P) 



(2) 



n{xyi<i>{^)\^ 



(n(x)) 



Where (. . .) is an ensemble average, 0(x) and Jn{x) are respectively: the phase 
and the absolute value of the wave function \i/(x) describing the whole system. As 
a coherence length Is of the system we take the full-width at half maximum of the 
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Figure 2. (Color online) Behavior of the density fluctuations at the center of the trap 
versus interaction for 1000 atoms. In (a) fluctuations divided by the density and in 
(b) divided by the square of density. 



gi{—x,x). Analogously the width of the atomic cloud is a full- width at half maximum 
of |\l'(x)p. In recent papers we showed that in both repulsive [6] and attractive [5] gas 
two regimes exist. One between the zero temperature and Tph in which the size of the 
condensate is smaller than the coherence length. The second one, above Tph, in which 
the opposite condition occurs - a quasicondensate regime. Such a situation is shown in 
Figure [H for weakly interacting repulsive gas. To check how important are fluctuations 
of the phase in the behavior of Z,^ we compared gi with: 

^i(-x,x) = (e-**(-^)e*'^(^)) (3) 

Function cji unlike gi accounts for phase fluctuations only. In Figure [1] we see that 
in our regime of parameters gi{—x,x) is nearly the same as gi{—x,x), so density 
fluctuations do not contribute to the first-order correlation function regardless the sign 
of the interaction. Moreover, density fluctuations are very similar for both attractive 
and repulsive cases. In Figure [2] we show S'^n/n (square of the density fluctuations 
divided by the density of atoms) at the center of the trap for changing g. We see 
that for the repulsive gas the fluctuations are smaller than for the attractive one. The 
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reason is that for the repulsive case the mean-field interparticle interaction energy is 
positive and proportional to the density squared. Because of this, the system prefers even 
distribution of the density rather than concentration of atoms. In the attractive gas we 
have the opposite situation. The interaction energy is lower for local bunching of atoms. 
Nevertheless, we don't see any qualitative difference between positive and negative g, 
as the dependence is nearly linear. What is more: the value of the fiuctuations divided 
by the square of the density is the same for every g what is presented in Figure [2j 
Finally, we see that coherence properties as well as density fiuctuations are very similar 
for attractive and repulsive gas. Therefore next, we take a closer look at the properties 
of the ideal gas. 

3. Ideal gas 

Detailed exact calculations are possible for the ideal gas. In the grand canonical 
ensemble results can be found in [T6l Sj |25] . 

We know that 0j(x) in ([8]) are the eigenstates of harmonic oscillator. In the 
canonical ensemble we used the recursion derived in [26] to calculate the partition 
function of A^ atoms: 

Zo =1, 

Zi(/3) =^exp(-/3e,), (4) 

V 

1 ^ 

where e^ is a single particle energy of the state i^, A^ is a number of atoms. Knowing 
Zat we can calculate probability of finding n atoms in a state of energy e^. Probability 
of finding at least n atoms there is: 

P3^(n|iV) = e--^-%^. (5) 

Therefore: 

PM^) = e-^^'-^^^^^ - e-(n+i)/3^.^^^^z!izi. (6) 

Zn Zn 

Then we can find the average number of atoms in the v mode: 

N, = Y.PMN)n. (7) 

n 

For deeper understanding of coherence properties lets go back to the definition ([2]) of 
gi{—x,x). We know that the one-body density matrix has the form: 

p(x, y) = {^*{x)^{y)) = E §'PK^)My), (8) 

where Ni/N and (/'i(x) are eigenvalues and eigenvectors of p{x,y). Hence, 

_ E.=of0:(-a:)0»(x) 
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Figure 3. (Color online) Comparison of ID and 3D ideal Bose gas in a harmonic trap. 
Populations of the ground state and the first excited state in both geometries for 1000 
atoms are shown. 



We see a strong connection between the coherence length of the system and the spectrum 
of one-body density matrix. Existence of a few comparable eigenvalues of the one-body 
density matrix signifies fragmented condensate [27]. It follows from the Penrose-Onsager 
criterion that the quasicondensate transition is caused by "sticking" of one-body density 
matrix eigenvalues - the situation when Nq is comparable to A'^i. 

There are two important remarks which enable us to gain better intuition before 
further analysis. First, when all atoms are in the condensate, this means that A^o = N, 
then gi{—x,x) is everywhere equal to one so we have a fully coherent system. On the 
other hand if all eigenvalues are equal then gi{—x, x) ~ 6{x — x') and coherence length is 
equal to zero. To show a notable difference in the form of density matrix spectrum in ID 
and 3D for the ideal gas we present in Figure |3] a population of two lowest modes (A^o, 
A'^i respectively) as a function of temperature. We see a qualitative difference between 
ID and 3D. In 3D Nq is much bigger than A^i all the way up to T^. Additionally the 
"sticking" of A^o and A'^i takes place for almost vanishing Nq, Nq/N <^ 1, i.e. almost at 
Tc. On the other hand, in ID we see the "sticking" for lower relative temperatures and, 
what is important, for much higher values of Nq/N. 

Using the intuition gained, we analyze more precisely one, two and three dimensions. 
The most intuitive way of comparing temperature effects for different regular geometries 
is to present temperature in Tc units: 

hoo N 



ID :T, 



2D ■.Tr= —- 



hu ( N V'^ 
C(2)j 



3D:r, 




(10) 



A^ 
03), 



1/3 



For 2D and 3D T^ is equal to the critical temperature in thermodynamical limit. In 
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Figure 4. (Color online) Ratio of populations of the ground state and the first excited 
state for ID, 2D and 3D in canonical ensemble (points) for N^/N ~ 0.2. Comparison 
with grand canonical results is presented (lines). Note a good agreement in 2D and 
3D. Grand canonical results for ID are calculated numerically using (fT2)) . In the inset 
is shown, that the discrepancy between both ensembles in ID exists even for huge 
number of atoms. 



ID the phase transition doesn't occur. This is why only characteristic temperature is 
used jl3]. However to make a comparison for different aspect ratios with finite number 
of particles it is more convenient to keep a relative population of the condensate Nq/N 
constant rather than the corresponding relative temperature. 

In Figure m we show Ni/Nq for different numbers of atoms in 1, 2 and 3 dimensions. 
All results are for Nq/N = 0.2. In ID 'sticking' is definitely much stronger than in 2D 
and 3D. Results in canonical ensemble are presented up to 1600 atoms - the number used 
in the experiment [3]. Because of numerical limitations the ratio of Ni/Nq for large A^ 
can be obtained in the grand canonical ensemble only. We see a good agreement between 
canonical and grand canonical ensembles in crossover range of A^ for 2D and 3D. 

The agreement is worse in ID. This is because the ensembles are equivalent in the 
thermodynamic limit only. Results of both ensembles (inset in Figure H]) approach each 
other on the logarithmic scale. Even for extremely large number of atoms the ratio 
Ni/Nq in ID is above 0.1, thus one would expect quasicondensation even for very large 
one dimensional system. Evidently, the finite size corrections in ID systems are quite 
important. The asymptotic behavior is not reached even for A^ ~ 10^^ particles. 

Moreover, the grand canonical approach allows also for finding the asymptotic value 
of the ratio of two dominant eigenvalues of the one-body density matrix. The quantum 
partition function for a system of noninteracting bosons in the grand canonical ensemble 
may be written as follows: 

Z{z^T) =Y.Y. ■■■X{z-'^M-Exnx/kBT), (11) 

no=Oni=0 A 

where Ex is the energy of the single particle state, and z = exp (fi/kBT) is a fugacity. For 
atoms trapped in a harmonic potential with frequency u we have in ID: E\ = huX, 2D: 
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Ex = huj{Xx+\y) and in 3D: Ex = huj{\x+\y+\z)- Knowing that (A^) = z^\nZ{z,T) 
we get an expression for the average number of atoms in the system: 

where we have extracted the term corresponding to the lowest energy Ex = 0. In 2D 
and 3D we can safely set 2; = 1 in the sum over all excited states (fT2|) . Next, we can 
expand the formula (fT2l) considering only leading terms in a small parameter (]hu |28] . 
it leads to the expression for the atom number: 

where D = 2 (3) in 2D (3D). Remembering that we are interested in results for the 

constant value No = CN = j^^we get z = j^ and ^ = ( ^^p^ )^^^- In the ID in 
the thermodynamic limit we have 

N = -^-Hl-z)^. (14) 

1 — z nuj 

Finally, we get that -^ is equal to ^^,^ ~ j. . The asymptotic behavior of Ni/Nq can be 
obtained from the relation: 

lim — - = r— r . (15) 

A^^oo No 1 - ze-^'^'^ z ^ ' 

Finally, in ID we get the following scaling of the 'sticking' ratio: 

lim — - ~ -7—r--T, (16) 

iv^ooiVo In(iV) ^ ' 

while in higher dimensions we have: 



that is 



and 



Ari/2 



for D 



for D 



iV2/3 

The sticking ratio has totally different behavior in the asymptotic limit depending on 
the dimensionality of the system. In higher dimensions it goes to zero with dimension- 
depend power while in ID the decay of the sticking ratio is logarithmically slow. This is 
the reason while in ID system there is always a finite range of temperatures below the 
quantum degeneracy temperature where occupation of the higher orbitals of the one- 
body density matrix is relatively high. The spectrum of one-body density matrix enables 
us to look at the dependence of T^^ (the temperature when the width of the system is 
equal to the coherence length) as a function of A^. In Figure [5t^a) for better comparison 
we divided T^h by Tc (ITOj) . Once more a significant difference between ID and higher 
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Figure 5. (Color online) Parameters of the system at the pomt of the transition from 
true condensate to quasicondensate versus number of atoms. In (a) temperature of the 
transition is presented, in (b) fraction of atoms in the ground state at the transition 
temperature. 



dimensions is seen. For 2D and 3D the temperature of quasicondensation is really near 
Tc whereas for ID the temperature Tph is much less than T^. Because of the ambiguity of 
the Tc definition in ID we looked at N^f^ (a number of atoms in the ground state at the 
temperature Tp^). In Figure [5] (b) we show N^y^/N as a function of A^. The ratio N^f^/N 
is a fast decreasing function tending to zero in 2D and 3D, while for ID it decreases 
much slower exceeding the value 0.5 in the whole presented region. Having analyzed 
ID, 2D and 3D symmetric cases we study the ratio of Ni/Nq and N2/NQ versus trap 
aspect ratio u}z/uj±. In Figure[6]we present results for fixed: A^ = 1000 and Nq/N = 0.4. 
This way we can analyze transition between different geometries of the system. When 
ojz <^ u}± the system is nearly ID, when uz ~ w± is nearly 3D, and when Uz ^ uj± is 
nearly 2D. Once more we see a clear difference between dimensions. However, we do 
not observe any sharp transition between them. The left vertical line corresponds to 
ksTph = huj± and the right to ksTph = huz, so to the situations when system can be 
considered respectively as one and two dimensional. We see that these points agree with 
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Figure 6. (Color online) Relative number of atoms in the ground state and in the first 
excited state versus the aspect ratio ojz/(^± for N = 1000 and Nq/N = 0.4. The left 
vertical line denotes point when ksTph = huj± and the right one when ksTph = TiiOz- 



the beginning of the flattening of Ni/Nq and N2/NQ. Moreover 2D case is much more 
similar to 3D than to ID. In experiments we never have exactly symmetric systems, 
however relatively broad range of aspect ratios around Uz = uj± corresponds to systems 
which can be considered as 3D because of the flattening of Ni/Nq in this region. We 
checked that this region gets broader with increasing A^. 

Summarizing, we have shown that quasicondensation phenomenon is strongly 
related to the dimensionality of the system rather then to interactions. Using the sticking 
ratio Ni/Nq as the quasicondensation criterion we found the asymptotic behavior of this 
quantity in the ideal gas case. The ID asymptotic of the sticking ratio differs drastically 
from asymptotic in higher dimensions. Vanishing of the second largest eigenvalue of the 
one-body density matrix is logarithmically slow in ID. Some decreasing of the correlation 
length in 2D and 3D can be also observed but this effect results from the finite number 
of atoms in the system and occurs just below the critical temperature. 
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